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Abstract. This paper is a report on a computer check of some important 
positivity properties of the Hecke algebra in type H4, including the non- 
negativity of the structure constants in the Kazhdan-Lusztig basis. This 
answers a long-standing question of Lusztig's. The same algorithm, carried 
out by hand, also allows us to deal with the case of dihedral Coxeter groups. 

1 Statement of the problem 

1.1. Let W be a Coxeter group, S its set of distinguished generators, and 
denote < the Bruhat ordering on W. Denote Ti. the Hecke algebra of W over 
the ring of Laurent polynomials A = Z[v, v ], where v is an indeterminate. 
We refer to [3] for basic results about Coxeter groups and Hecke algebras; 
we just recall here that Ti is a free A-module with basis {t y ) y& w, where the 
algebra structure is the unique one which satisfies 



tg.ty 



t sy if sy > y 

(v - v~ l )t y + t sy \isy<y 
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(here we use t y = v~ l ( y 'T y , where the T y satisfy the more familiar relation 
T s .T y = (q — l)T y + qT sy when sy < y, with q = v 2 .) Then there is a unique 
ring involution i on H such that i(v) = v" 1 , and i(t y ) = t~_ x . 

The Kazhdan-Lusztig basis of H is the unique family (c y ) ye w such that 
(a) i(c y ) = c y and (b) c y = t y + ^2x<yPx,ytx, with p x<y G v^Zf-i; -1 ]; in 
particular we find that c s = t s + u ~ 1 for all s G S. It turns out that 
Px,y = o ! 'w _! Wp !t « is a polynomial in g, the Kazhdan-Lusztig polynomial 
for the pair x, y. For any pair (x, y) of elements in W, write 

C-x-C-y — ^ ^ h x y z C z 

(in other words, the h xy z G A are the structure constants of the Hecke 
algebra in the Kazhdan-Lusztig basis). 

A number of the deeper results in the theory of Hecke algebras depend 
on positivity properties of the polynomials P XtV and h x ^ y ^ z . More precisely, 
consider the following properties : 

P±: the polynomials P x<y have non-negative coefficents; 

P2: the P xy are decreasing for fixed y, in the sense that if x < z < y in 
W, the polynomial P x>y — P ZtV has non-negative coefficents; 

P3: the polynomials h x>y z have non- negative coefficients. 

Properties P\ and P3 are basic tools in the study of Kazhdan-Lusztig 
cells and the asymptotic Hecke algebra; they have been proved in (Hj and 
[5] for crystallographic W using deep properties of intersection cohomology 
(see the remarks in section 3 of Lusztig jj] for the case where W is infinite.) 
Property P2 is proved for finite Weyl groups in jlj, using the description of 
Kazhdan-Lusztig polynomials in terms of nitrations of Verma modules. 

None of these geometric or representation-theoretic interpretations are 
available in the non-crystallographic case, and the validity in general of the 
positivity properties above are among the main open problems in the theory 
of Coxeter groups. Let us concentrate on the case where W is finite. It is 
easy to see that for the validity of P\— P3 we may reduce to the case where 
W is irreducible. 

1.2. The case where W is dihedral is simple enough to be carried out by 
hand; we have included the computation of the h x>y z in section^ This leaves 
us only with the two groups H3 and H4, and since the former is contained 
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in the latter, the only case we need to consider is the Coxeter group of type 
H4, of order 14400, with Coxeter diagram 

12 3 4 
o o o o 

5 

Despite its rather modest size, the group H4 poses a redoubtable compu- 
tational challenge, even for present-day computers. It is all the more re- 
markable that property P\ was checked already in 1987 by Dean Alvis by 
explicitly computing all the Kazhdan-Lusztig polynomials. Quite a feat with 
the hardware of the time! Unfortunately, Alvis's programs have never been 
made available; to my knowledge, the only publicly available computer pro- 
gram capable of carrying out this computation is my own program Coxeter 
0, which does it in about one minute on a modern-day personal computer. 
This still leaves open properties P2 and P3; the main purpose of this paper 
is to report on a computation, carried out as one of the first applications of 
version 3 of Coxeter, by which we prove : 

1.3 Theorem. — Properties P2 and P3 hold for the Coxeter group of 
type H4. 

1.4 Corollary. — The fifteen conjectures labelled PI P15 in Chapter 
14 of Lusztig JEy hold for the Hecke algebra of the Coxeter group of type 
H4, and in fact for the equal parameter Hecke algebra of any finite Coxeter 
group. 

Proof. - Lusztig shows in Chapter 15 of 8^ that in the equal parameter 
case (which is automatic in type H4), the fifteen conjectures all follow from 
Pi and P3 . In view of our earlier remarks on the reduction to the irreducible 
case, and of the dihedral computation in section |1J the same argument can 
be applied to any finite Coxeter group. 

2 Methodology 

2.1. The verification of P\ and Pi is straightforward : one simply runs 
through the computation of the P x>y for all x < y in the group. In fact, 
from the well-known property P x>y = Psx,y whenever sx > x, sy < y (cf. jS] 
(2.3.g)), and the analogous property on the right, for Pi it suffices to consider 
the cases where LR(x) D LR(y), where we denote L(x) = {s E S \ sx < x} 
(resp. R(x) = {s G S \ xs < x}) the left (resp. right) descent set of x, and 
LR(x) = L(x)\J R(x) C SJJS"; we call such pairs (x,y) extremal pairs. 
Taking also into account the symmetry P X}V = P x -i „-i, there are 2 348 942 
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cases to consider, which are easily tabulated by the program. The non- 
negativity of the polynomials is checked as they are found. For P2, one also 
easily reduces to extremal pairs (x,y) and (z,y). The tough computation 
is for P3; here there are a priori 14400 3 = 2 985 984 000 000 (almost three 
trillion!) polynomials h XiViZ to be computed, and the only obvious symmetry 
is h XtV:Z = hy-i x -i z -i (but, as explained below, we don't even use that.) 

2.2. The algorithm used in the computation is simple. For a fixed y, we 
compute the various c x .c y by induction on the length of x, starting with 
c e .c y = Cy, where e denotes the identity element of W. To carry out the 
induction, we choose any generator s £ S such that sx < x, and write : 

c x = c s .c sx - ^2 fi(z,sx)c z (1) 

z<sx 
sz<z 

where as usual /j,(x, y) denotes the coefficient of v^ 1 in p X)V (which is also 
the coefficient of degree \(l(y) — l(x) — 1) in P x , y , and in particular is zero 
when the length difference of x and y is even.) 

Then we may assume that c sx .c y is already known, and similarly for the 
various c z .c y , so we are reduced to multiplications of the form c s .c u , for s £ S 
and u £ W. When su > u this is read off from formula (1) with x = su; 
and when su < u one simply has c s .c u = (v + v ~ 1 )c u (see for instance [3]). 

The information that is required for this computation is encoded in the 
VF-graph of the group; once this graph is known, everything else just involves 
elementary operations on polynomials. Actually, it is obvious that the h Xt y, z 
are in fact polynomials in (u+t> -1 ), so they are determined by their positive- 
degree part; it is this part that we compute and keep in memory. The only 
complication arises from the fact that we need to be careful about integer 
overflow; it turns out in fact that for H4 all the coefficients of the h XjVjZ fit 
into a 32-bit unsigned (and even signed) integer, but not by all that much: 
the largest coefficient that occurs is 710 904 968, which is only about six 
times smaller than 2 32 = 4 294 967 296. 

2.3. From the computational standpoint, this procedure has a number of 
desirable features. First and foremost, once the VF-graph of the group has 
been determined, the problem splits up into 14400 independent computa- 
tions, one for each y, so we can forget about the computation for a given 
y when passing to the next (this would not be true if we tried to use the 
symmetry h x ^ V)Z = h y -\ x -i g -i .) This is advantageous in terms of memory 
usage, could be used to parallelize the computation if necessary (it turns 
out that it hasn't been), and also means that the computation can be harm- 
lessly interrupted (either voluntarily or involuntarily) , at least if its progress 
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is recorded somewhere : basically, the only penalty to pay for picking up 
an interrupted computation is the recomputation of the lU-graph, which 
takes about half a minute. This is very valuable for computations running 
over several days, where there is always the risk of a system crash or power 
failure. 

On the other hand, it is essential that for a fixed y, the full table of h x ^^ z 
be stored in memory. In practice, there are many repetitions among these 
polynomials; so we store them in the form of a table of 14400 2 = 207 360 000 
pointers. Initially, this is the main memory requirement of the program; it 
is interesting to note that the cost doubles, from about 800 MB to about 1.6 
GB, when we pass from a 32-bit to a 64-bit computer. It turns out that for 
a fixed y, the additional memory required to store the actual polynomials is 
small, and never exceeds 300 MB. So the full computation runs confortably 
in 2 GB of memory, and barely exceeds 1 GB on a 32-bit machine. 

It is much more difficult to try to write down a full table of all the 
polynomials that occur as h X;y>z . I have done this a number of times, but 
with the memory available on the machines to which I have had access, it 
has been necessary to split up the computation in about a hundred pieces 
to avoid memory overflow in the polynomial store, to keep the correspond- 
ing files in compressed form to avoid overflowing the hard disk, and then 
to merge those one hundred compressed files into a single compressed file, 
eliminating repetitions. At some point I have needed to store about 30 GB 
of compressed files — not something administrators are very happy about! 
In view of these difficulties, I have chosen not to make that version of the 
program available for the time being. 

2.4. The computation has been done a number of times (writing files of 
all the distinct polynomials): at the Ecole Polytechnique, Centre de Calcul 
Medicis, Laboratoire STIX, FRE 2341 CNRS, on several computers, includ- 
ing a Compaq Alpha EV68 and an AMD Opteron server with 4 and 8 GB 
of main memory respectively, where it has taken about 5 days of CPU time, 
and twice at the Universite Lyon I, Maply, UMR 5585 CNRS, on a Xeon 
processor with 4 GB of main memory, where it took 80 hours (not counting 
the final file merging pass, which took another 80 hours or so.) The pro- 
gram as presented here has run on our 2.7 GHz AMD Opteron server at the 
Institut Camille Jordan, using less than 2 GB of memory, in a little under 
85 hours. 

On the technological side, it seems that the time was just right for this 
computation : it makes full use of the 3Ghz processors, at least 2GB central 
memories, and 100+ GB hard disks that are found on typical low-end servers 
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today. It would still be beyond the grasp of most present-day personal 
computers, however, although that, too, is changing fast! 

3 Verification 

3.1. Let's come now to the thorny issue of verification : what is the amount 
of trust that can be put in a result like this ? An obvious prerequisite is the 
availablility of the source code of the program that carries out the computa- 
tion; this may be downloaded at http : / / igd.univ-lyonl . f r/~ducloux/ coxeter/coxeter3/positiv: 

This is actually the source code of an especially modified version of 
Coxeter3. All the extra code is contained in the file special. cpp; all the 
other files are identical to the ones in Coxeter3. 

3.2. In addition to those already available in Coxeter3, the following com- 
mands are defined : 

- klplist : prints out a list of all the distinct Kazhdan-Lusztig polyno- 
mials which occur in the group (so that in particular, one may re-check 
property P\, although this was done already by the computer in the 
course of the computation); 

- decrklpol : checks property P2; 

- positivity : checks the non-negativity of the h XiV z ; 

- cycltable : prints out a table of the c x .c y for a fixed y; 

- cprod : prints out a single product c x -c y ; 

For type H4, on a decent server, the cycltable command should not 
take more than two minutes; cprod should usually take less than a minute 
(Note that the first call will take longer than subsequent ones, because the 
VF-graph must be computed the first time.) So these commands give local 
access to the multiplication table of the Hecke algebra, thus opening up the 
"black box" a little. 

The positivity command for type H<± can also be executed through 
the little stand-alone command coxbatch that I have included; this will 
run the computation in the background. It writes any errors in error_log, 
and records the progress of the computation in positivity_log. After a 
successful run, error_log should be empty, and positivity_log should end 
with the line : 

14399: maxcoeff = 710904968 

(elements of the group are numbered from to 14399.) 
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3.3. As was explained in 12,21 the computation is entirely elementary once 
the VK-graph of the group is known. The trust that one may place in this 
ingredient should be rather high, in my opinion, as it is computed with an 
algorithm which in simpler cases has been checked against other programs, 
and which even for H4 has been checked against other algorithms for the 
same computation. (The latter check is perhaps more convincing than the 
former, for it may very well happen that some of the nastiest configurations 
occur for the first time in type H4, or even only in type H4, as it is such an 
exceptional group.) 

For the actual non-negativity check, we are as far as I know in entirely 
uncharted territory. However, the code for the computation of the h Xt y, z 
from the VF-graph is really quite simple, so it can rather convincingly be 
checked by inspection. Another check is as follows : the order in which the 
computations are performed depends on the choice of a descent generator 
for x. By default, we always choose the first such generator in the internal 
numbering of Coxeter; however, it is easy enough to replace this by other 
"descent strategies" (for instance, choosing the first generator in some other 
ordering.) This will lead to a very different flow of recursion. I have done 
this replacing the first descent generator by the last one, and obtained the 
same output file, which is as it should be. 

4 The case of dihedral groups 

4.1. Let us now show how the algorithm described in 12.21 can be carried 
out by hand in the case where W is dihedral. This is not difficult and may 
well be known to experts, but I haven't been able to find the results in the 
literature. In Chapter 17 of [B], Lusztig computes the h X y >z for an infinite 
dihedral group with unequal parameters. I have not been able to determine 
if his formulas can be specialized to yield the statement in Proposition 14.41 
(although of course our computations for the infinte case are very similar to 
(and much simpler than) those in |Sj.) 

4.2. Assume first that W is infinite dihedral. Let S = {s%, S2}, and for each 
i > denote 



where in each case there are i terms in the product. When we need indeter- 
minate generators, we will also use notation such as [s,t,i> for s ^ t in S. 



[l,2,i> = sis 2 • • • 
<i, 1,2] = ...s 1 s 2 



[2, l,i> = s 2 si • • • 
<i,2,l] = . ..s 2 si 
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Denote for simplicity c\ = c si , c 2 = c S2 . It is well-known that for dihedral 
groups all the Kazhdan-Lusztig polynomials P x>y are equal to 1 for x < y; 
incidentally, this proves that Pi and P 2 are trivially verified. It follows that 
formula [2.21 (1) reduces to 

C\.C2 = c Sl . S2 = C[ 12)2 > 

ci-C[ 2) i,t> = C[i,2,i+l> + c [i,2,i-l> for % > 1 

and of course ci.C[i )2i i> = (i> + u" 1 ) C[ 12j j> for i > 0, and similar formulas 
for the left multiplication by c 2 . 

4.3. Now fix an element of the Kazhdan-Lusztig basis, that without loss 
of generality we may assume to be of the form c < fc i 2 ], with k > 0. Let s 
be the first generator in <l,2,k], so that c^i^i = cr Sj4) ^ ;> , and let £ 7^ s 
be the other element of S. We wish to compute the various c<£ s ,t]' c <k,l,2]> 
c <j,t,s]- c <fc,i,2]i as * varies. Note that since the c^i^i, 3 > 0, form a basis 
of a left cell representation of Ti, it is a priori clear that these products will 
be ^4-linear combinations of the c <Jjl 2 ]. Start with the c <i s ft- c <k,l,ih wnere 
we may assume i > 0. For the first two values of i we get 

c <l,s,i]- c <fc,l,2] = c t-C[ Sy t,k> = c [t,s,k-l> + c [t,s,k+l> = c <fc-l,l,2] + c <fc+l,l,2] 
c <2,s,t]- c <fc,l,2] = C s- C t- C <A;,1,2] = c <fc-2,l,2] + ^ c <k,l,2] + c <fc+2,l,2] 

when k > 2, and 

{ c t- c <i,i,2] = C <2,1,2] , , 

for k = 1 
Cs.Ct.02. = C<1,1, 2 ] + C< 3 ,1,2] 

( Ci - C<2 ' 1 ' 2l=C<1 ' 2 ' 11+C<3 ' 2 ' 11 fork = 2 

{ C s -Ct-C<2,l,2] = 2c <2,2,l] + C <4,2,1] 

Now applying the procedure from 12.21 yields the recursion formula: 

c <i,s,t]- c <k,l,2] = Cr-C<i-l, s ,t]' c <k,l,2] ~ c <i-2,s,t]- c <fc,l,2] for i > 2 

where r € {s, t} is the first term in c < j s i ]. It follows easily that the non-zero 
terms in c <i s ^.c <k i 2 ] all correspond to indices j of the same parity, which 
changes when we go from one i to the next. Consequently, all these terms 
have a first generator not equal to r. If we write the coefficients in rows, we 
see that the coefficient at position j in row i > 2 is obtained by adding the 
coefficients at positions j — 1 and j + 1 in row i — 1, and subtracting the 
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coefficient at position j in row i — 2 . For example, when k = 3, the table 
looks like this (with dots indicating zeroes) : 

k-2 fe-1 k k + 1 k + 2 k + 3 k + 4 k + 5 

i=l • 1-1 

i = 2 1 2 1 

i = 3 • 2-2 ■ 1 

i = 4 1 2 2 • 1 

j = 5 • 1-2 ■ 2 ■ 1 

Due to the exceptional multiplication rules ci.c 2 = c <2j 2,i] an d c 2 .ci = 
c <2,i,2]; the first 1 that would appear in column should be omitted (this 
occurs for i = k); another way of stating this is that we should run the 
algorithm in the half-plane j > 0, and omit undefined terms. In this form, 
the procedure is valid for all values of k > 0. Note that for i > k the rows 
have constant length k + 1, the next one being just the previous one shifted 
one unit to the right. 

The case of the c <i t s yc <k i 2 ] is similar but simpler; it appears as Propo- 
sition 7.7 (a) in Chapter 7 of |Hj. Here all the non-zero coefficients are equal 
to v + v _1 , and we have a second order recursion as previously. Again the 
coefficients start expanding pyramidally, there is a "reflection" at column 0, 
and then we get a band of constant width k. 

The final result may be stated as follows: 

4.4 Proposition. — Let W be an infinite dihedral group, and adopt the 
notation from \4.'A Let k > 0, and let s be the Brst generator in c < fc i 2] (i.e., 
s = 1 if k is even, s = 2 if k is odd.) Let t be the other element in S. Then 
we have : 

c <i,s,t]- c <k,l,2] = 

c <k-i,l,2] + 2c <fc-i+2,l,2] + • • • + 2 C<fc+i-2,l,2] + c <fc+i,l,2] (0 <i < k) 

2 C<2,1,2] + • • • + 2c <2 fc_2,l,2] + C< 2 fe,l,2] (« = k ) 

C<i-fc,l,2] + 2 C<j-jfc+2,l,2] + • • • + 2c <i+fc-2,l,2] + C<i + fc : l )2 ] {i > k) 

(where for k = 1 the entry for i = k should be interpreted as c < 2 i i i 2] ), &nd 
similarly 

c <i,t,s]- c <k,l,2] = 

(v + v' 1 ) c <fe _j +lj i )2 ] + ... + (v + v' 1 ) c <fc+ i_i )lj2 ] (0 < i <= k) 
(v + v- 1 ) c <i _ fc+l!l)2] +... + (« + v- 1 ) c <i+fc _ 1)lj2 ] (i >= k) 
with of course similar formulas when c <fe l 2 ] is replaced by c <k 2 ,i] 
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4.5. Now consider the case where W is a finite dihedral group of order 2m, 
m > 2. Then of course c < j jl 2 ] is defined only for % < m, and moreover the 
action of a and c 2 on c <m 12 ] = c <mj2) i] is given by 

Cl-C< m ,l,2] = C 2 .C <m) i )2 ] = (W + W _1 )c <mjli 2] 

The consequence is that the recursive pattern described above gets modified 
starting from i = m + 1 — k. For that row, the expression c <m _i j i j 2] + 
c <m+i,i,2] that would have been obtained by applying the appropriate c s to 
c <m,i,2] should be replaced by (t> + v~ 1 )c <m ^^]- The net effect is that the 
algorithm splits in two independent parts: one one hand we run the same 
algorithm as for the infinite case, but this time within the strip < j < m; 
on the other hand, starting from i = m — k, we add a term of the form 
a-ic <m ,i,2], with a m - k = 1, a rn - k+1 = v + v' 1 , a m - k+2 = v 2 + 2 + v~ 2 , and 

at = (v + v _1 )ai_i - ai_ 2 = w d + 2w d " 2 + . . . + 2v~ d+2 + 

for d = k + i — m > 2. The procedure goes on until i = m, at which 
point the row from the first algorithm has disappeared altogether, and only 
the multiple of c <m l 2 ] is left. For example, when m = 9, k = 6, and the 
c <i,i,2]- c <fc,i,2] ( we have s = si in this case), we get the following table for 
the first part of the algorithm : 

fc-5 fc-4 fc-3 fc-2 fc-1 fc fe+1 k + 2 




and row 9 is zero. The obvious symmetry in the shape is a general fact. 

Note that if we look at the Z-basis of H afforded by the v d c w , d G Z, and 
give v d c < j i 2] degree j + d, then the sum of the coefficients of c <iyS ^ .c^ 12 i 
in each degree is the same for the finite and the infinite cases; whatever 
is missing from the infinite picture is exactly reflected in the coefficient of 
c <m,i,2]- So a concise statement of the result is as follows : 
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4.6 Proposition. — Let W be dihedral of order 2m, m > 2. Then the 
formula? from Pro position \4.4\ rem am valid, except that one must have i < m, 
and that any expression of the form 

c <m+d,l,2] + c <m-d,l,2] d> 

that can be taken out of the formula should be replaced by (v d +v~ d )c <m i^] ■ 

4.7 Example. — Pursuing the earlier example where m = 9 and k = 6, 
and taking for instance the case where i = 6, putting together the sixth row 
in the above table and the expression for the coefficient a^, which corresponds 
to d = 3, we get: 

C <6,1,2]- C <6,1,2] = 

2C< 2 ,1,2] + 2 C<4,1,2] + C< 6 ,1,2] + + + 2v^ + V~ 3 ) C <9 ^ 2] (1) 

The corresponding expression for the infinite group would be 

2c <2,l,2] + 2c <4,l,2] + 2 C<6,1,2] + 2c <8i i i2 ] + 2c < i 0) l,2] + C <12,1,2] 

from which we recover (1) by taking out out one copy of c <6j1j2 ] + c <12j i i 2] 
and two copies of c <81)2 ] + c<io,i,2]- 

It is easy to get many other examples from the program — of course 
playing with the program is how the above statements were found in the 
first place. 

5 Questions 

5.1. On looking at the lists of polynomials which occur as h Xt y jZ , one imme- 
diately notices that they are not only non-negative, but have a much stronger 
positivity property : if we denote d the degree of h Xj y :Z , then v d h X: y tZ is a 
polynomial in q = v 2 , which is unimodal (recall that this means that the 
coefficients increase up to a point, which in this case has to be the middle be- 
cause of the symmetry of the h X) y^ z , and decrease from there.) In the course 
of the computation, the program checks unimodality for all h x<y ^ z , and prints 
an error message on error_log in case of failure. Hence unimodality holds 
for the Hecke algebra of type H4. 
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5.2. For Weyl groups, there is one case where it is easy to prove that 
unimodality holds : viz. the case where y = u>o is the longest element in the 
group. From the properties of the <lr preorder it is clear that A.c wo is a 
two-sided ideal in H; so we may write c x .c Wo = h x c Wo , where h x = h X)VjQ)Wo . 
Now it is clear that t s .c WQ = v.c wo for all s, hence t x .c m = v l ( x 'c Wo , and 

z<x 

from which it follows immediately, using the expression of the intersection 
homology Poincare polynomial in terms of Kazhdan-Lusztig polynomials (jSJ 
Theorem 4.3.) that v l<yX ^h x is equal to the Poincare polynomial of the in- 
tersection (hyper)cohomology of the Schubert variety X x . The unimodality 
then follows from the so-called hard Lefschetz theorem. As far as I know, 
the unimodality property for general h x> y tZ is an open question, even for 
Weyl groups. 

5.3. Clearly, all the results about the Hecke algebra of type H4 which are 
stated in this paper point to the fact that there is a hidden geometry here 
that is begging to be discovered. Hopefully, the facts about this geometry 
which the program opens up will help us understand what is going on, and 
serve as a guide towards the solution. I should be very happy if this turns 
out to be the case. 
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